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Abstract 

The order submission and cancelation processes are two crucial aspects in the price formation of stocks traded in 
order-driven markets. We investigate the dynamics of order cancelation by studying the statistical properties of inter- 
cancelation durations defined as the waiting times between consecutive order cancelations of 22 liquid stocks traded on 
the Shenzhen Stock Exchange of China in year 2003. Three types of cancelations are considered including cancelation 
of any limit orders, of buy limit orders and of sell limit orders. We find that the distributions of the inter-cancelation 
durations of individual stocks can be well modeled by Weibulls for each type of cancelation and the distributions 
of rescaled durations of each type of cancelations exhibit a scaling behavior for different stocks. Complex intra- 
day patterns are also unveiled in the inter-cancelation durations. The detrended fluctuation analysis (DFA) and the 
multifractal DFA show that the inter-cancelation durations possess long-term memory and multifractal nature, which 
are not influenced by the intraday patterns. No clear crossover phenomenon is observed in the detrended fluctuation 
functions with respect to the time scale. These findings indicate that the cancelation of limit orders is a non-Poisson 
process, which has potential worth in the construction of order-driven market models. 
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1. Introduction 



Order submission and cancelation are two central processes in the price formation of stocks traded in order-driven 
markets. Understanding their statistical regularities are crucial in the study of stock market microstructure theory 
and the construction of order-driven models lUSSl- For limit orders, there are three attributes: order direction (or 
order sign indicating buy/sell), order price, and order size. The statistical properties of these quantities in the order 
submission process have been extensively studied including the long memory of order signs |lll|4|], the distribution of 
relative prices d S 0, 0, S H US U 1 12 1 . the long memory of relative prices known as the"diagonal effect" fljs , 3 1 



22, 



and the power- law distribution and long memory of order sizes and trading volumes 

n Q s m uti m in, E 

iL 13211 . In contrast, the statistical regularities of the order cancelation process are less 
studied. When constructing the empirical behavioral model of order-driven stock markets, the conditional probability 
of order cancelation on three factors has been determined Jl]]. In this work, we attempt to understand the order 
cancelation dynamics by investigating the inter-cancelation durations that are the waiting time between consecutive 
cancelations. 
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There are numerous studies on the waiting time distributions of diverse financial quantities. Generally speaking, 
the waiting time is defined as the time interval between two successive financial events. When the financial events 
are defined, one is able to determine the series of waiting times. One importance topic is about the return intervals 
(or recurrence intervals). On the one hand, the recurrence intervals between financial volatilities exceeding a certain 
thresho ld q have been carefully studied, and numerous phenomena have been unveiled 1 33, 34, 35, 3^ 37, 38, 3^ 40i 
4ll l42l I43U44II . On the other hand, econophysicists have devoted to the study of recurrence intervals between large 



positive or negative price returns, which has important implications on risk estimation |45, 4^ 47, 48, 49 

Another important topic is about the intertrade durations, which are defined as the waiting times between con- 
secutive transactions of an equity. The importance of this topic is due to the fact that intertrade durations contain 
information contents of trading activity and have crucial relevance to the microstructure theory ifsoi 51 , 52]. In main- 



stream finance, the celebrated autoregressive conditional duration (ACD) model was proposed to model the intertrade 
durations with temporal correlation and other financial variables ||53, 54]. Alternatively, in the econophysics com- 
munity, the continuous-time random walk (CTRW) formalism has been adopted to deal with the intertrade durations 
and price dynamics ]|55, 5^ 57, 58, 5^. Empirical investigations of the intertrade durations from different equity 
markets report that the probability distributions might be described by power laws IHol 61], modified power laws 



I67[] . In addition, a scaling behavior 
which is however less conclusive 



157, 59I1. stretched exponentials (or WeibuUs) 1 25, 62, 63, 64, 65, 6^ 67], stretched exponentials followed by power 
laws implying that the transaction process is non-Poisson. However, statistical tests reject the hypothesis that 

the intertrade durations are distributed according to an exponenti al JtoI ItHi or a power law [ |66l l. Hence, the Weibull 
seems to be the very form of intertrade duration distribution 1621 163[ I64l l65i l6i 
can be observed in the distributions of rescaled intertrade durations ll64[ l66l l67l 

lEl. 

The cancelation of limit orders has certain impact on the price formation. Especially, when the orders at the best 
bid or ask prices on the order book are canceled, the bid-ask spread widens and the mid-price changes. In some order- 
driven models, the cancelation is assumed to follow a Poisson process ll72ll . Although this is a good approximation, we 
find in this work that the inter-cancelation durations are distributed according to a Weibull rather than an exponential 
and long-term correlated. These findings show that the cancelation of limit orders is a non-Poisson process, which is 
useful in constructing more realistic order-driven market models. 

The rest of this paper is organized as follows. In Section |2] we briefly describe the data sets investigated and the 
basic statistics of limit order cancelations. Section |3] investigates the empirical distributions of the inter-cancelation 
duration. In section H] the temporal correlations and the multifractal nature of the inter-cancelation durations are 
studied based on the (multifractal) detrended fluctuation analysis. Section[5]summarizes and concludes. 



2. Data description 

Shenzhen Stock Exchange (SZSE) is an order-driven market. There were three different periods on each trading 
day in the SZSE before July 1, 2006, namely, the opening call action (9:15 am to 9:25 am), the cooling period (9:25 am 
to 9:30 am), and the continuous double auction (9:30 am to 1 1:30 am and 13:00 pm to 15:00 pm). More information 



of interest about this market can be found in the literature 111 1LI73L 17411 . This study is based on the ultrahigh-frequency 



data of 22 liquid stocks traded on the SZSE, which contain all order placements and cancelations in the whole year 
2003. Since part of the data of stock 000002, 000027, 000063, 000088, 000089 are missing, the data from 39-th to 
167-th trading days of these stocks are discarded. We only take into consideration the cancelations occurring in the 
continuous double auction, and the interval across the pausing period is excluded. Since the orders not executed during 
the continuous double auction period will be canceled automatically by the system, no inter-cancelation duration will 
be calculated overnight either 

Inter-cancelation duration is defined as the waiting time between two consecutive cancelations (in units of second). 
The resolution of time is 0.01 second. When the (/ -1- l)-th cancelation order arrives, the value t, = f,+i - f, is recorded 
as the i-th inter-cancelation duration, where f, is the time when the /-th cancelation occurs. The cancelations in the 
market can be classified into two types based on their order directions, namely the cancelations of buy orders (buyer- 
initialed cancelations) and the cancelations of sell orders (seller-initialed cancelations). By analogy with the definition 
of T,, the duration (r?) between successive cancelations of buy (sell) orders is defined as 



2 



(1) 



where (f?) is the time stamp when the i-th buyer-initialed (seller-initialed) cancelation occurs. There are less than 
2% cancelations occurring at the same time for each stock. 

Table[T]lists the number of cancelations (A^, N^,N'^), the number of simultaneously happened cancelations (A^o, N^, A^^) 
and the average inter-cancelation duration ((t), (t''}, (t")) for the cancelation of all orders, buy orders, and sell orders 
for the 22 Chinese stocks under investigation. We note that there are several trivial relations in this table stating that 

(1) the total number of cancelation is the sum of the cancelation numbers of buy orders and sell orders 

N ^N'' + N-\ (2a) 

(2) there are less simultaneous cancelations of distinct buy/sell Umit orders than all limit orders 

A?o>A?o+^o. (2b) 

and (3) the mean duration of consecutive cancelations of buy orders or sell orders is greater than the mean duration of 
all cancelations 

(T^) > (T> and <T^) > (T>. (2c) 

These relations hold for all stocks. 

Table 1: Descriptive statistics of the 22 Chinese stocks studied over the whole year 2003. Since part of the data of stock 000002, 000027, 000063, 
000088, 000089 are missing, the data from 39-th to 167-th trading days of these stocks are not taken into consideration. N is the number of 
cancelations, Nq is the number of simultaneous cancelations, and (t> is the average inter-cancelation duration for the 22 Chinese stocks. 



Code 


All cancelation 


Buyer-initiated cancelation 


Seller-initiated cancelation 


A^ 


M) 




N" 












000001 


591944 


11178 


5.77 


317015 


5141 


10.75 


274929 


4012 


12.38 


000002 


78378 


1220 


20.05 


34577 


726 


45.14 


43801 


391 


35.66 


000009 


371225 


3440 


9.27 


183804 


1183 


18.68 


187421 


1347 


18.30 


000012 


221023 


1853 


15.37 


1 14662 


659 


29.40 


106361 


849 


31.82 


000016 


119408 


1069 


28.58 


60219 


444 


56.34 


59189 


525 


57.29 


000021 


310027 


3193 


11.05 


157174 


1633 


21.75 


152853 


1068 


22.32 


000024 


88187 


1068 


38.11 


42593 


741 


78.17 


45594 


268 


73.01 


000027 


71705 


1006 


21.74 


33058 


545 


46.77 


38647 


392 


40.09 


000063 


60496 


1136 


26.05 


25681 


660 


60.68 


34815 


414 


44.98 


000066 


216984 


1554 


15.77 


110289 


643 


30.90 


106695 


577 


31.82 


000088 


13778 


625 


112.29 


6861 


410 


220.82 


6917 


212 


219.50 


000089 


43893 


1129 


35.30 


20909 


486 


73.45 


22984 


605 


66.98 


000429 


73173 


527 


46.15 


36999 


311 


90.50 


36174 


164 


92.15 


000488 


66439 


1424 


51.02 


32585 


583 


102.72 


33854 


789 


99.24 


000539 


54037 


3790 


61.95 


26950 


3000 


120.76 


27087 


748 


122.47 


000541 


39562 


580 


85.82 


19715 


218 


169.87 


19847 


351 


167.93 


000550 


252471 


4696 


13.57 


122865 


3236 


27.75 


129606 


1042 


26.25 


000581 


56726 


1798 


59.56 


27236 


1284 


122.47 


29490 


483 


112.81 


000625 


255080 


5399 


13.23 


123361 


3263 


27.24 


131719 


1753 


25.49 


000709 


129721 


1267 


26.17 


65704 


602 


51.20 


64017 


517 


52.61 


000720 


30767 


3805 


103.55 


16558 


2490 


175.61 


14209 


1309 


207.04 


000778 


90936 


1022 


37.38 


43576 


532 


77.06 


47360 


426 


71.36 
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3. Empirical probability distributions 



3.1. Empirical distributions 

In this section, we study the empirical probabiHty distributions of the inter-cancelation durations for individual 
stocks. The three types of the empirical distributions of the inter-cancelation durations t,t* and t' for the 22 individual 
stocks are plotted in Fig.[Tla-c). A very important feature is that all these distributions are not exponential. In addition, 
it is not evident that these distributions have power-law tails by eye-balling. A horizontal comparison of the three plots 
reveals that the the three distributions P{t), Pii^), and P(t^) are very similar 
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Figure 1: (Color online) Empirical probability distributions of the inter-cancelation durations for all limit orders, buy limit orders and sell limit 
orders, (a-c) The empirical probability density of inter-cancelation durations for all cancelations, buyer-initialed cancelations and seller-initialed 
cancelations for the 22 stocks, (d-f) ScaHng in distributions of the rescaled inter-cancelation durations, (g-i) Fitting the probability density of the 
rescaled inter-cancelation durations of the ensemble of all the 22 stocks to the WeibuU and the ^-exponential distributions using the maximum 
Ukelihood estimatoion and the nonhnear least-squares regression. 

Similar to the situation of intertrade durations 1164 , 16711 . we conjecture that the inter-cancelation duration 
distributions for individual stocks are Weibulls, 



PwBL(T)=y8ff-V-le-'^ 



(3) 



Following Refs. Il66ll67ll . we also try to model the distributions using ^-exponentials 



l+(l-q){—) 



(4) 



We have fitted each curve in Fig.fTla-c) using the above two equations based on the maximum Ukelihood estimation 
(MLE) and the nonlinear least-squares estimation (NLSE). The estimated parameters for the curves in Fig. [Ha) are 
depicted in Tabled The rm.s. of the fitting residuals ;i'wBL of the Weibull distribution is less than;^'qE of the q- 
exponential distribution for 18 stocks when we use the maximum likelihood estimation. When we use the nonlinear 
least-squares regression, only one stock has x^'Rh > XqE- We thus conclude that the inter-cancelation durations are 
better modeled by the Weibull distribution. For the waiting times between successive cancelations of buy orders and 
sell orders, we obtain very similar results and the Weibull outperforms the ^-exponential as well. 



Table 2: Estimated parameters of the Weibull distribution (3) and the (^-exponential distribution ^4) based on the maximum likelihood estimation 
(MLE) and the nonlinear least-squares estimation (NLSE), respectively. The columns ^wBL i^^'iXqE 'he rm.s. values of fit residuals. 



Code 






MLE 










NLSE 






a 


P 




q 


XqE 


a 


P 


;tWBL M 


q 


XqE 



000001 0.41 0.67 0.42 0.24 1.67 0.51 0.18 

000002 0.47 0.71 0.15 0.30 1.57 0.30 0.34 
000009 0.42 0.69 0.42 0.26 1.62 0.54 0.31 
000012 0.35 0.66 0.31 0.19 1.71 0.25 0.22 
000016 0.48 0.71 0.17 0.32 1.54 0.27 0.32 
000021 0.39 0.68 0.29 0.23 1.65 0.39 0.18 
000024 0.44 0.69 0.16 0.28 1.60 0.33 0.29 
000027 0.42 0.69 0.24 0.26 1.61 0.30 0.27 
000063 0.36 0.64 0.39 0.18 1.81 0.42 0.26 
000066 0.38 0.67 0.29 0.22 1.67 0.34 0.24 

000088 0.56 0.72 0.40 0.41 1.45 0.62 0.44 

000089 0.51 0.73 0.15 0.35 1.49 0.24 0.40 
000429 0.42 0.68 0.19 0.25 1.66 0.31 0.31 
000488 0.41 0.67 0.43 0.24 1.68 0.55 0.28 
000539 0.39 0.63 0.29 0.21 1.79 0.43 0.35 
000541 0.51 0.71 0.14 0.34 1.52 0.36 0.36 
000550 0.33 0.64 0.34 0.17 1.78 0.30 0.19 
000581 0.44 0.69 0.23 0.27 1.62 0.39 0.34 
000625 0.34 0.65 0.34 0.18 1.73 0.24 0.15 
000709 0.39 0.67 0.22 0.23 1.66 0.20 0.26 
000720 0.23 0.52 0.26 0.08 2.30 1.06 0.17 
000778 0.45 0.70 0.18 0.29 1.59 0.25 0.33 



0.47 0.85 0.49 1.26 1.14 

0.57 0.60 0.50 1.27 0.77 

0.55 0.54 0.53 1.24 1.14 

0.49 0.66 0.43 1.30 1.14 

0.56 0.74 0.52 1.25 0.71 

0.47 0.81 0.46 1.28 1.05 

0.54 0.62 0.48 1.28 0.86 

0.53 0.76 0.46 1.29 0.82 

0.52 0.37 0.44 1.31 1.34 

0.50 0.71 0.41 1.33 0.96 

0.63 0.51 0.55 1.25 0.79 

0.61 0.59 0.55 1.24 0.63 

0.55 0.49 0.48 1.28 0.96 

0.53 0.47 0.46 1.30 1.19 

0.58 0.30 0.51 1.27 1.20 

0.58 0.60 0.54 1.24 0.73 

0.47 0.58 0.42 1.31 1.29 

0.57 0.46 0.51 1.26 0.94 

0.44 0.90 0.40 1.33 1.08 

0.52 0.67 0.47 1.28 0.98 

0.45 0.43 0.30 1.45 1.13 

0.56 0.62 0.55 1.23 0.86 



Table |2]illustrates that the estimated values of each parameter for different stocks are close to the mean, especially 
for p. In addition, the shapes of the empirical distributions shown in Fig. [Ua-c) are very similar to those of the 



intertrade durations 1164 166ll67ll . To the best of our knowledge, a scaling behavior in the rescaled intertrade duration 



distributions was first reported by Ch. Ivano v et al. t641 with further evidence provided by Jiang et al. 116711 . In order to 
check if the distributions of inter-cancelation durations also have a scaling behavior, we apply the following rescaling 
scheme to each stock 

T t/ct, P o-P, (5) 

where cr is the standard deviation of the durations of a given stock. Figure[Tld-f) plots the rescaled probability erf (r) 
as a function of the rescaled duration r/cr for the three types of inter-cancelation durations. We find in each plot that 
all the 22 curves collapse onto a single curve, which implies a scaling form 

ctP{t) = p(t/ct), (6) 

where p is the scaling function. 

Since the rescaled distributions of inter-cancelation durations for different stocks show a nice scaling, it enables 



5 



us to treat all the rescaled inter-cancelation durations from different stocks as an ensemble to gain better statistics. We 
use g, and g'^ to denote the variables t/ct, jcr^ and r^/cr* from 22 different stocks, respectively. Figure [Ttg-i) 
illustrates the three empirical probability densities of the three types of rescaled intercancelation durations. Again, we 
use the Weibull and the ^-exponential to model each distribution pir/cr) based on the maximum likelihood estimation 
and the nonlinear least-squares estimation. The fitted curves are also drawn in Fig.[Tlg-i), and the estimated parameters 
are listed in Table |3] We find that;t'wBL < XqE for all the cases except for the inter-cancelation durations of sell limit 
orders g^ based on the maximum likelihood estimation where ;ifwBL is slightly greater than;t'qE. Therefore, the Weibull 
distribution is a better model for the three scaling functions. We also observe that the three p values are close to each 
other for either the maximum likelihood estimation or the nonlinear least-squares regression. It is interesting to note 
that the yS values are also close to those for the intertrade durations ll67ll . 



Table 3: Estimated values of parameters (a,j3, q, fi) by means of MLE and NLSE. ;^'wbl and^qE stand for the rm.s. of fit residuals. 









MLE 










NLSE 








a 


P 


^WBL A* 


q 




a 




/t'WBL A* 


q 




g 


0.40 


0.67 


0.21 0.23 


1.67 


0.30 


0.21 


0.49 


0.70 0.55 


1.23 


1.13 




0.38 


0.65 


0.27 0.21 


1.75 


0.43 


0.23 


0.50 


0.45 0.51 


1.25 


1.28 


g' 


0.34 


0.64 


0.26 0.18 


1.78 


0.24 


0.11 


0.42 


0.84 0.47 


1.27 


1.26 



We now focus on the curves in Fig. [TJg) fitted using the Weibull distribution. It is found that the maximum 
likelihood estimation curve fits the bulk of the distribution very well in the range [10^^ 5], which accounts for 99.1% 
of the sample. However, it is unable to capture the large durations. In contrast, the nonlinear least-squares regression 
curve fits the tail quite well but deviates markedly the bulk of the distributions. 

3.2. Conditional distributions of inter-cancelation durations 

We now investigate the conditional distribution of rescaled inter-cancelation durations on the value of its preceding 
duration. All the rescaled durations for different stocks constitute an ensemble set Q. We sort the set Q in an increasing 
order and partition it into eight non-overlapping groups of the same size: 

8 

e = U 2" (7) 

where Q, f] Qj - (p when / j, and gi < gj when gt e Qi, gj e Qj, i < j. We then estimate the empirical conditional 
probability distribution P(g\go e Qi) = P{g(t)\g(t -1)6 Qi), which is the probability density of the rescaled inter- 
cancelation durations g(t) whose preceding value go - g(t - 1) belongs to Qi. 

The eight empirical conditional PDFs are plotted in Fig. |2a). Assuming that ; > j, we note that P{g\gQ e Qi) < 
P{g\go e Qj) for small g and P{g\gQ e Qi) > Piglgo 6 Qj) for small g. In other words, large durations tend to follow 
large durations and small durations are prone to follow small durations. Figure |2b) shows the dependence of the 
conditional mean duration (glgo) with respect to go. It is shown that the conditional mean duration increases with go, 
which is consistent with the outcome of Fig. |2 a). This phenomenon indicates that there is short-term memory in the 
inter-cancelation durations. 

We recall that the conditional distributions of the intertrade durations with respect to different preceding durations 
almost collapse onto a single curve and the conditional mean of intertrade durations does not changes with the pre- 
ceding durations fE^. Therefore, although both the intertrade durations and the inter-cancelation durations exhibit 
scaling in the distributions and the scaling functions are both Weibulls with very close exponents, the difference in the 
behavior of conditional durations unveils that transactions and cancelations follow different dynamic processes. 

4. Memory effects 

4.1. Intraday pattern 

Intraday patterns exist in many high-frequency financial variables in the empirical studies. It is necessary to 
investigate the intraday patterns in the inter-cancelation durations for each stock. To obtain the intraday patterns of 

6 




Figure 2: (Color online.) (a) Conditional probability density P(g\go € Qi). (b) Conditional mean inter-cancelation durations (g\go) with respect to 
go = <g : g 6 Qi)- 



inter-cancelation durations, we segment the continuous double auction of each trading day into 240 successive 1-min 
intervals. For a given stock, the inter-cancelation durations are then averaged within each trading minute to create a 
minute-by-minute series as follows. 



1 ' 



(8) 



where Njj represents the number of inter-cancelation durations in the j-th interval in the i-th trading day, ta is the 
inter-cancelation duration of an order which is canceled in the y-th minuite, and T,y is the average duration of the y-th 
interval in the i-th trading day. The mean inter-cancelation duration at a given minute of a trading day is calculated as 
follows. 



(9) 



where is the number of trading days. Figure |3]plots {T}j as a function of the intraday time for four randomly chosen 
stocks (000001, 000027, 000581, 000709). Complex intraday patterns are observed, which share some analogues with 
those of intertrade durations 117 511 . 
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Figure 3: (Color online) Intraday pattern of inter-cancelation durations for four different stocks traded on the SZSE duiing the calendar year 2003. 
The dots show the cross sectional mean values of inter-cancelation durations over the whole trading year at that minute. The continuous curves are 
the polynomial fits to the data. From the curves, a crude inverse f-shaped pattern can be observed. 



Since the inter-cancelation duration shows an explicit intraday pattern, we define the socalled adjusted inter- 
cancelation durations by removing the intraday patterns from the original data, 

f, = t,/{t)j. (10) 

We will analyze both the original data and the adjusted data in the following, which allows us to determine if the 
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intraday patterns have influence on the temporal correlations. 



4.2. Long-range dependence 

To detect the long-range dependence of the inter-cancelation durations, the detrended fluctuation analysis (DFA) 
is utilized, which is able to extract long-range power-law correlation in non-stationary time series ll76ll77ll . The DFA 
procedure consists of following steps. For a given inter-cancelation duration series {t,|! - 1,2,..., A^}, the cumulative 
summation series y, should be first calculated as follows. 



1,2,..., A?. 



(11) 



Then we use A^, disjoint intervals with the same size s to cover the series y. Since the length of the series need 
not be a multiple of the size of the interval s, the whole series y, may not be completely covered by A^, intervals. For 
compensating the remain part, we can use another A^, intervals to cover the series from the other end of the series. In 
each interval, a polynomial is used to calculate the local trend function y by least-squares regressions. In this section, 
linear functions are used in the fitting procedure. The local detrended fluctuation function rk{s) in the ^-th interval is 
defined as the r.m.s. of the fitting residuals: 



(12) 



where h is the k-\h interval. Thus the overall detrended fluctuation is then estimated as follows 

2N, 



(13) 



To determine the scaling behavior of the fluctuation function, we vary the scale s in the range of [20, A^/4] (scale s > 
N/4- is excluded since the detrended fluctuation F2{s) becomes statistically unreliable). Thus a power-law relationship 
between the overall fluctuation function F2(s) and the interval size s can be expected as follows. 



F2(s) ~ s". 



(14) 



where H stands for the DFA scaling exponent. Practically, we can plot F2{s) as a function of s on double logarithmic 
scales to measure // by a linear fit. 
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Figure 4: (Color online) Log-log plots of the overall fluctuation function Fi{s) with respect to the interval size i for four different stocks. The 
open circles and squares stand for the original and adjusted data, respectively. The solid lines are the power-law fits to the data. The curves for the 
adjusted data have been shifted downwards for clarity. 

Figure |4] shows the log-log plots of the overall fluctuations ^2(5) as a function of the interval size s for the 
four randomly selected stocks. Excellent power-law dependence is observed for each curve with the scaling range 
spanning more than two orders of magnitude. No clear crossover is observed in the present case, which is however 
quite common for other financial quantities |78, 79, S^l- We also find that the two curves in each plot parallel to each 
other, which means that the intraday patterns have minor influence on the long-term power-law correlations. 
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In Table ID the DFA scaling exponents H for all the 22 stocks are listed. All the DFA scaling exponents are 
significantly greater than 0.5, indicating that the inter-cancelation durations have strong long memory. The DFA 
scaling exponents H of both the original data and the adjusted data do not differ much, which confirms that the intraday 
patterns have little influence on the long-range dependence of inter-cancelation durations. This finding provides 
further evidence that the cancelation process is a non-Poisson process. 



Table 4: Chai'acteristic parameters of the long-term memory and the multifractal nature of the inter-cancelation durations. The first column gives 
the stock codes. The second column lists the mean number of cancelations {Nj) in one trading day. The third and the fifth columns list the DFA 
scaling exponents H\ of the original data and H2 of the adjusted data. The rest two columns list the width of multifractal spectrum Aa of the 
inter-cancelation durations for the 22 stocks. 



Code 


{Nt) 




Orig 


inal data 




Adjusted data 




H 


1 










H2 


Aa 


000001 


2476 





91 


+ 





02 


0.76 





94 


+ 0.00 


0.88 


000002 


327 





93 


+ 





02 


0.76 





96 


+ 0.00 


0.89 


000009 


1546 





97 


+ 





05 


0.53 


1 


01 


±0.00 


0.49 


000012 


928 





83 


+ 





03 


0.58 





86 


+ 0.00 


0.59 


000016 


497 





89 


+ 





04 


0.67 





92 


+ 0.00 


0.60 


000021 


1291 





90 


+ 





02 


0.74 





93 


±0.00 


0.69 


000024 


370 





90 


+ 





03 


0.54 





93 


±0.00 


0.46 


000027 


300 





90 


+ 





02 


0.72 





94 


±0.00 


0.71 


000063 


254 





89 


+ 





05 


0.74 





91 


±0.00 


0.70 


000066 


903 





89 


+ 





04 


0.59 





92 


±0.00 


0.53 


000088 


57 





80 


+ 





17 


0.61 





83 


±0.00 


0.54 


000089 


185 





92 


+ 





04 


0.47 





95 


±0.00 


0.39 


000429 


305 





85 


+ 





03 


0.82 





88 


±0.00 


0.69 


000488 


277 





88 


+ 





07 


0.59 





91 


±0.00 


0.43 


000539 


226 





86 


+ 





09 


0.81 





88 


±0.00 


0.76 


000541 


164 





85 


± 





06 


0.60 





89 


±0.00 


0.57 


000550 


1051 





87 


+ 





05 


0.73 





90 


±0.00 


0.68 


000581 


236 





84 


+ 





05 


0.89 





87 


±0.00 


0.76 


000625 


1071 





89 


+ 





05 


0.75 





91 


±0.00 


0.65 


000709 


542 





92 


+ 





05 


0.69 





94 


±0.00 


0.71 


000720 


129 





84 


+ 





20 


0.74 





85 


±0.00 


0.69 


000778 


379 





92 


+ 





05 


0.64 





95 


±0.00 


0.55 



4.3. Multifractal nature 

We now investigate the multifractal nature of inter-cancelation durations by using the multifractal detrended fluc- 
tuation analysis (MF-DFA) 18 ill . The first several steps of the MF-DFA procedure are essentially the same as the DFA 
procedure, and the multifractal detrended fluctuation is defined as the following form, 

C , 2iV, Y 
W^l^^gMin > (15) 

where the moment order q varies in real number except for q -Q in which Fq{s) is defined as 

Fo{s) = exp 1^ g Inms)]^ . (16) 
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For each given value q, we vary the value of s in the interval [20, as in DFA. If the series is long-range power 
correlated, a power-law relation is expected between the detrended fluctuation function Fq{s) and the size s. 



F,{s) ~ s''^^\ (17) 

where h(q) is the MF-DFA scaling exponent. Note that when q -2, h{2) is nothing but the DFA scaling exponent H. 
The singularity strength a and its spectrum /(a) can be obtained based on the Legendre transform, 

a = qh'{q) + h(q) 
f(a) = q^h'(q)+l ' 

To investigate the multifractality of the inter-cancelation durations, the fluctuation functions Fg(s) are calculated 
for each stock. Since the size of each time series is finite, the estimate of Fg{s) will fluctuate remarkably for large 
values of \q\, especially when the scale s is large. We focus on ^ € [-6, 6] to obtain reasonable statistics in the 
estimation of Fg(s). We show the detrended fluctuation functions of the original data and the adjusted data for the 
four typical stocks. Again, nice power laws are observed for all the curves. 




Figure 5: (Color online) Dependence of the detrended fluctuation functions F^(.v) of the original data (upper panel) and the adjusted data (lower 
panel) with respect to the scale s for the four randomly chosen stocks. The solid lines are the best power-law fits to the data in the con'esponding 
scaling ranges. 



For each given q, the MF-DFA scaling exponent h(q) can be estimated by the linear regression between ln[F^(s)] 
and In s. The resultant h(q) functions are shown in Fig.|6ta), which decrease with the increase of q. Figure|6lb) illus- 
trates the multifractal spectra for the four stocks. It is evident that the inter-cancelation durations possess multifractal 
nature, on which the intraday patterns have minor impacts. 

The strength of the multifractal nature can be quantified by the width of the singularity spectrum 

Aa = ttmax - ffmin- (19) 



A large value of Aa corresponds to stronger multifractality. The values of the singularity width Aff for the 22 stocks 
are calculated and listed in Table |4] These Aff values are significantly different from zero as a hallmark of strong 
multifractality. 



5. Conclusion 

We have investigated the inter-cancelation durations calculated from the limit order book data of 22 liquid Chinese 
stocks traded on the SZSE in the whole year 2003. It is found that the probability densities decrease with the increasing 
of the inter-cancelation durations. The 22 empirical densities of the rescaled durations collapse onto a single curve, 
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Figure 6: (Color online) Multifractal analysis of the inter-cancelation durations for four typical stocks. Panel (a) shows the MF-DFA scaling 
exponent h(q) with respect to q, and panel (b) shows the corresponding multifractal spectra f(a). 

showing a nice scaling pattern. This scaling behavior is also observed in the distributions of buyer-initiated inter- 
cancelation durations and seller-initiated inter-cancelation durations, which implies that there are common features in 
the cancelation behavior of market participants. 

We then model the three ensemble densities of rescaled inter-cancelation durations for all cancelations, buyer- 
initialed cancelations and seller-initialed cancelations by using the Weibull and the Tsallis ^-exponential. By using 
the maximum likelihood estimation and nonlinear least-squares estimation methods, it is found that all the three kinds 
of inter-cancelation durations can be well modeled by the Weibull function. We also study the conditional distribution 
of rescaled inter-cancelation durations which follows a certain set of rescaled inter-cancelation durations. It is found 
that large (resp. small) durations tend to follow large (resp.) durations. It means that the inter-cancelations exhibit 
short-term memory. 

In addition, the intraday pattern, the long memory, and the multifractal nature of the durations are also investigated. 
An inverse U -shaped intraday pattern is revealed for each stock. The DFA and the MF-DFA methods are applied to 
both the original data and the adjusted data after removing the intraday patterns. We show that the intraday patterns 
have minor influence on these correlation properties. The results show that the durations possess both the long memory 
and the multifractal nature. 
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